function List = subnormalFPN(beta,p,L,U)
pList = [];
e = L;
l = beta^(e-p+1);
while (l < beta^L)
    pList = [pList,l];
    l = l + beta^(e-p+1);
end
List = [-fliplr(pList),0,pList];